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Pulsars: What we see 

radio: 
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• Pulse peaks are narrow 
Negligible energy budget 


PWNe feed by dense plasma 
• Energy goes there 
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Pulsar Magnetosphere: Emission regions 

"Friction" in the machine 
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Pulsar Magnetosphere: Large scale view 
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Pulsar Magnetosphere: Outer parts - Wind Zone 

Numerical simulations of Force-Free magnetosphere 
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(Kalapotharakos et al. 2012) 

very good match with analytical solution of Bogovalov {1 999) 
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Pulsar Magnetosphere: Structure of inner parts 

Numerical simulations of Force-Free magnetosphere 



(Spitkovsky 2006) 
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Pulsar Magnetosphere: The global circiut 
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Current density in the polar cap 

current and charge densities are not refated: j = ncjC 


from 2D modeling of an aligned rotator 
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Current density in the polar cap 

current and charge densities are not related: j ^jtj= ti G jC 

from 3D modeling of an inclined rotator 



[thanks to Xue-Ning Bai (Bai & Spitkovsky 2010)] 


L. Mestel, Yu. Lyubarsky, some others: “told you so!” 
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Physical Model: Plasma flow in an electric circuit 
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Boundary Conditions: Models for NS surface 


• [RS] 

particles cannot leave NS surface 
Ruderman & Sutherland (1975) model 

self-consistent model: Timokhin 2009, 201 0 


• [SCLF] 

particles freely leave NS surface Space Charge Limited Flow 
Arons & Scharlemann (1979) model 
(also Muslimov & Tsygan 1992, Harding & Muslimov 1998, 
and others) 

self-consistent model: Timokhin & Arons 2012 
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Numerical Model: Modeling from first principles 

Numerical Code: 

Particle acceleration Electric field Particle-ln-Ce.il 
Particles — > Photons -» Particles( Pairs) tonte Carlo 

Physical Ingredients: 


• 1 D Electrodynamics 

starting from 

0E ll „ , 

— _47rt^-n=,l 

• y-ray production : curvature radiation 

• pair creation : single photon absorption in magnetic field 
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Numerical code for cascade modeling 


■ | Enforce Gauss Lav | 
Adjust number of particles 
~ ~ ~ T ~ 

Apply particles 
homirUry conditions 
: i ; 


Inject pairs 


for each particle: 


if emits: 

for each photon: 

simple photons 

if absorbed: 


create a pair 


. »orte Carlo 
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Adjust number of pairs 
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Save Results 
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• Move particles . 

| 


■ Scatter p and j 

■Solve Maxwell’s Equations I 

Q t:^t + M 

• Collect fluxes 
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Modeling from the first principles: 
Particle acceleration Electric field 
Particles -» Photons -> Particles(Pairs) 
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RS: No particles supplied by the NS j = j C) 

Limit cycle: series of discharges 


particles ’ momenta p = -7 
c 
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Magnetosphere ► 


0.0 0.2 0.4 0.6 0.8 1.0 


X distance from NS 


• electrons positrons • y-rays 
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RS: No particles supplied by the NS j j G , 

Single discharge 


particles’ momenta p = -7 
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SCLF: Low energetic flow 0 < j/j C] < 1 

Formation of the low energetic flow 

V 

particles : momenta p= 7 
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Ad - Debye length x [Ad] dictance fivm NS 
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| — NS Magnetosphere ► 

0 10 20 30 40 50 


Dashed red line - oscillating solution in one-fluid approximation 
Mestel & Pryce (1985), Shibata (1997), Beloborodov (2008) 
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SCLF: Low energetic flow 0 < j/j G} < 1 

Current and charge density adjustment: beam and cloud 



particles* momenta p = -y 
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SCLF: Low energetic flow 0 < j/j C] < 1 

Maximum particle momentum as a function of £, = //jbj 
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SCLF: Discharges in super-GJ flow j/j C) > 1 


particles’ momenta p~-y 
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j — NS Magnetosphere 
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electrons positrons • y-rays * protons 
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Discharges in the return current j/j Q < 0 

The same for both RS and SCLF 
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particles* momenta p = — 7 
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j‘ outer" magnetosphere 

| no pair formation 
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Particle acceleration in SCLF regime 
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Particles' 

momenta: 

p~[0,10] 

electrons 


Polar cap 
currents 
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Particles' 

momenta: 

p— [0,10 7 ] 
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Particle acceleration in RS regime 


"hot" flow everywhere, 
pair production: 


near NS 0 < j/j Q , < 1 



Particles 1 

momenta: 

p~[0,10 7 l 
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photons 


Polar cap 
currents 
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Waves during discharge 

It did not escape our attention. . . 


RS__wave_xp_e_jp05 .mp4 
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Conclusions 


• Current density determines the plasma flow regime 

• Cascades are non-stationary. ALWAYS. 

• All flow regimes look different: mutiple components (?) 

• Return current regions should have particle accelerating zones in 
the outer magnetosphere: y -ray pulsars (?) 

• Plasma oscillations in discharges: direct radio emission (?) 
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